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Abstract 

We have obtained the solutions of linearized Shrddinger equation for spher- 
ically and axially symmetrical electrons density oscillations in plasma in the 
approximation of the self- consistent field. It was shown that in the center or 
on the axis of symmetry of such a system the static density of electrons can 
enhance, which leads to the increasing of density and pressure of ion gas. We 
suggest that this mechanism could be realized in nature as rare phenomenon 
called the 'fireball' and could be used in carrying out the research concerning 
controlled fusion. 

If a volume charge appears in electroneutral plasma, i.e. electrons density 
increases or decreases in some finite area, then, after an external influence 
is over, the oscillating process consisting in periodical changes of the sign of 
the considered volume charge is known to appear. It is necessary to remind 
that the frequency of this process ui p called plasma frequency is related to 
free electrons density in the plasma no by the formula 

2 4iie 2 n 



where e and m are the charge and the mass of the electron. 

We will assume electrons in plasma to be a quantum many body system. 
Such an assumption is due to the fact that, as it will be shown below, the 
solutions obtained have characteristic sizes of atomic order. It is known 
that the classical dynamics of iV interacting particles can be represented 
by the system of differential equations of motion in configuration space of 
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3iV-dimensions, or in 6iV-dimensional phase space, as well as by the system 
of partial differential equations in three-dimensional physical space and as 
the dynamics of singular material fields (see Refs. [1-3]). The transfer to 
the quantum-mechanical description is realized by changing the dynamic 
functions to Hermitian operators. In this case the dimensionality of the 
configuration space is conserved and the state of the system is completely 
defined by the wave functions in the 3iV-dimensional space. 

Shrodinger equation as it was shown in Refs. [3,4], could be presented 
in the form of the system of equations in physical space and having the 
same form as the equations of hydrodynamics. In these works the quantum- 
mechanical system of N particles with arbitrary masses and charges, inter- 
acting between itself by Coulomb forces and with external classical electro- 
magnetic field characterized by vector A and scalar ip potentials has been 
considered. 

In this approach the complex \1/ function introduced in three-dimensional 
space has the following form 



*(f f ,t) = yJn e (r,t)e*' T ™, (2) 

where |\1/| 2 = n e (f,t) is the density of electrons and \1/ function satisfies the 
partial differential equation: 

<9\I> 

iR - = ™. (3, 

The Hamiltonian in Eq. (3) is expressed in the following way 

1 ( h \ 2 ? 

H = - — ( -V — -A (f, t) ) + ecp + e 2 [ dh'G (r - 7) \^\ 2 + 6 (r,t), (4) 

where G (r — r f ) = pc^j. 

In the Eq. (4) the first two terms are the components of one electron 
Hamiltonian in the external electromagnetic field, the third term represents 
the potential energy of the electron in the self-consistent electrostatic field 
created by the whole electron system, with density of number of particles 
being equal to |^| 2 . The function 9(f,t) describing exchange interactions 
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between electrons has the following form 



0/ ^ A h 2 Afr f dp e (r,t) 2 r 3 , /■ ^ g 2 (r, r*, t) 



ro ro 



where p e (r, i) is the pressure of electron gas, g 2 (r, r*, t) is the correlation 
function. 

Therefore, in order to resolve exactly the considered problem one should 
take into account all terms in the Eq. (5). However, we will assume that 
the contribution of exchange interactions to the dynamics of free electrons in 
plasma is much smaller than that of the potential of self-consistent electro- 
static field. Then, in describing the dynamics of electron gas we will neglect 
the function 9 (r , t) . As it will be seen from further speculations, this rough 
approximation allows us to get some consequences which are close to those 
obtained from the consideration of the similar problem within the classical 
approach. 

Let us consider the neutral plasma, formed by the singly ionized gas, with 
the energy of electrons being more than the potential of ionization of this gas. 
We will suppose that this plasma possesses a spherical symmetry for density 
and velocities distribution of the electron and the ion gases. Taking into 
account the small mobility of heavy ions in gas compared to that of electrons 
we will suppose that the density of ions is the constant value rij (r,t) = Uq. 
We will also consider that in our case there are no external electromagnetic 
fields except those of positively charged ions. It is worth mentioning that the 
magnetic field in the system is equal to zero. The potential of self-consistent 
field created by the density of electron gas is represented by the formula 
(using spherical coordinates): 

U e = e "J dhGif - r)\^\ 2 =47re J ^ j x 2 \^(x, t)\ 2 dx. (6) 

ro r 

Similarly, for the potential ip of singly ionized gas with density of ions n^r, t) 
one has 

7dR f 

ip = -47re J —J x 2 rii(x,t)dx. (7) 
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Thus, taking into account Eq. (6) and Eq. (7), the Eq. (3) can be represented 
in the following way 

ihj- + — A^-47re 2 ^F(|^| 2 -n ) = 0, (8) 

where A = + 2 J^ is the Laplas operator in the spherical coordinate 

oo R 

system, F(. . .) = J §Jx 2 (. . .)dx. 

r 

Moreover, we demand the system to be electroneutral as a whole, i.e. the 
condition must be satisfied: 

1 f 2|,T,/_. ,m2j„ ™o 



lim / x 2 |^(x,t)|^x = ^. (9) 



We will search for a solution of the Eq. (8) in the form: 

$ = $ + xe -*"* . (10) 

Supposing that |\I>o| 2 = n and \x\ < |* |, we get |^ | 2 = ^o + ^o/ + |x| 2 , 
where / = xe~ lw * + x* e%ult - Taking into account that \x\ <C |^ |, we obtain 
|^o| 2 ~n + ^ /. 

Then, we substitute Eq. (10) and approximate expression for l^ol 2 i n 
Eq. (8). Having considered the complex-conjugate equation together with 
the obtained one, it was easy to get: 

huf + - 4vre 2 M/o (2* + /) F (/) = 0. (11) 

For the total linearization, it is necessary to suppose that 2^0 + / — 
2^o i n the third term of the equation involved. Then, we represent the 
function x through its real and imaginary parts: x — Xi + Hence 
/ = xi cos cut + x 2 sinu;t, and Eq. (8) can be divided into two independent 
similar equations for xi and xi'- 

h 2 

hooxn + ^ A Xn - 8vre 2 n F ( Xn ) = 0, n = 1, 2. (12) 
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One can find out the functions Xn = B n smir are the solutions of of the 
Eq. (12) if 7 satisfies the following dispersion relation 



urn 



1/2- 

1± 1 




(13) 



where u p is determined in Eq. (1). The positive part of the dependence of 
7 on the frequency cu is presented in the Fig. 1. It is worth noticing that if 
the values of B n are limited, the expressions for Xn satisfies the condition of 
system electroneutrality (9). 

From the expression (13) one can see that the frequency u = 2oo p is the 
critical value, since for frequencies less than 2u p , 7 becomes a complex value 
and under these circumstances the oscillations of the electron gas are absent. 
It is necessary to remind that in using the classical approach to the similar 
problem one gets the value u = u p for the critical frequency. 

For example, for the density tiq = 2.7- 10 l9 cm~ 3 , i.e. for completely singly 
ionized gas under the atmospheric pressure and when u = 2u p , the frequency 
of electron oscillations is: 



2 \ !/ 2 
7T \ Tim 



,= ^ = 2 ^1 ^9-10 13 ^. 



This frequency corresponds to the electromagnetic radiation in the in- 



frared range with the wavelength A = - ~ 3 • 10 cm. In this case 71 
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_l_ ^ 2m ^p j ' ~ -j-2.3 • 10 7 cm 1 , the size of central region 5, where the most 
intensive oscillations of the electron gas are observed, is equal to n/j m 
1.4 • 10~ 7 cm. 

It is necessary to remind that the exact expression for the density of 
electron gas in searching for the solution in the form of the Eq. (10) is pre- 

I 1 2 I 1 2 

sented as j^ol = ^0 + ^0/ + • Let us define in this expression fy f and 

1 1 2 

n e — \x\ + n o as the dynamic and static components of density respectively. 
In deriving the approximate linearized equation (12) the function \x\ 2 was 
supposed to be small and thus neglected. This procedure is not correct be- 
cause the integral i 7 " (|x| 2 ) for the function x = B^ 11 - is divergent. However, 
we assumed that the density of ion gas was constant throughout the volume 
because of the small mobility of heavy ions. This assumption is valid for 
the frequently oscillating dynamic component, but \x\ does not depend on 



1/2 
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71,2 = ± i^jr) i an( i f° r the lower branch 7 3i4 = ± 




time. Therefore, it is naturally to expect that under some conditions the 
negative volume charge, described by this function, will be compensated (or 
neutralized) by removing positive ions, that will result in local changing of 
density and pressure of the ion gas. 

Indeed, it can be shown that for our case the condition of static stability 
of the ion gas is (see, for instance, Refs. [5,6]): 



where k is Boltsman constant, T is the ion gas temperature, is the ions 
density, n e is the static electrons density. 

We demand that the ions density is equal to static electrons density with 
high level of accuracy rii ~ n e = no + B 2 sm r2 7r ■ Then, from Eq. (14) we get 




(14) 



o 
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the following inequality: 4^^- 757 <C 1. Having substituted the expression 
for 7ij in the last formula, we obtained the condition of the neutralization: 
k ™2n ^ 1- Taking into account that ions density in the center of the system 
should be equal to n c = B 2r y 2 , this condition can be rewritten in the following 

way: 

For instance, for T = 1000°-K" and the value 7 = 2.3 • 10' 1 cmT x which was 
obtained above we get n c <C 500no- In deriving the solution of our problem 
we assumed that n c <C no- Thus, we can conclude that, if the condition (15) 
is satisfied, the ion charge is undoubtedly compensated by the supplementary 
static electron charge and the divergence in the integral F(|^| 2 — rii) can be 
eliminated. 

Hence, the approximate, linearized theory trends to describe the pressure 
enhancement of the ion gas in the center of symmetrically oscillating electron 
gas. 

It is worth mentioning that along with spherically symmetrical solution 
of the Eq. (12), there is at least one axially-symmetrical one which has the 
form: 

= nl 12 +5J ( 7 r)e-^, 

where Jo is the zero-order Bessel function. In this case the dispersion relation 
takes the same form as Eq. (13). All consequences obtained for spherically- 
symmetrical oscillations are valid and for this case too. 

While considering nonlinear Shrodinger equation (8), it can be seen that 
along with the components of frequency u, the terms which do not depend 
on time as well as on frequencies 2tv, 3u and etc. appear. One can make 

sure of this representing the solution of the Eq. (8) in the form: \&(r, t) = 

00 

^0 + E *k(r,t), where ^ = B^'^^^. 
k=i r 
In Ref. [7] we put forward the hypothesis that in this case the static 

density of electron gas could attain great values. Density and pressure of ion 

gas are also great in the center of such systems. Thus, along with the heating 

of the central area resulted from intensive movement of electrons, there will 

be a possible running of nuclear fusion reactions of corresponding nuclei in 

gas. 

This process seems to support the existence of the enigmatic natural 
phenomena, known as 'fireball'. If vapors of water which are present in 
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atmosphere contain deuterium in the amount of m 1/5000 under normal 
conditions, the running the nuclear fusion reactions will release the energy, 
which support the oscillations of electron gas and prevent the recombination 
of plasma. Taking into account small sizes of the central (active) regions and 
small amount of deuterium in atmosphere, it is possible to use the term 'mi- 
crodose' thermonuclear reaction for process in question. Axially-symmetrical 
oscillation of the electron gas is likely to appear as very seldom observed type 
of a 'fireball' in the form of shining, sometimes closed cord [8]. Uncomplicated 
calculation shows that energy released in deuterium nuclei fusion in ldm 3 of 
vapors of water (the average size of a 'fireball') has the value of about 1MJ, 
that corresponds to evaluations of energy released by some observed 'fireballs' 
(see, for example, Refs. [8,9]). 

Groups and separate researchers developing the problem of the controlled 
fusion are suggested to pay attention to self-consistent, radially and axially 
oscillating electron plasmoids as a base models to self-supported thermonu- 
clear reactions. The authors of the present article have certain experience 
in generating ball-like plasma structures which proves the correctness of the 
chosen model. 
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